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NOTATION 



English Letters 



A 


Angstrom unit 


a 


equilibrium speed of sound 


b a 


partition function of argon atom 


V 


partition function of argon ion 


C 


speed of light 


D 


diameter of cylindrical leading edge of model (l/2 inch) 


e i 


specific ionization energy 


e t 


internal energy per unit mass due to translation 


h 


specific enthalpy 


h P 


Planck' s constant 




Gladstone-Dale constant 


k 


Boltzmann' s constant 


L 


test section width 


M 


Mach number 


M 

s 


incident shock Mach number 


m A 


mass of argon atom 


m 

A + 


mass of argon ion 


ra 

e 


mass of electron 


n a 


number density of atoms 


“a* 


number density of ions 


N 

e 


number density of electrons 
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NOTATION (Continued) 



n index of refraction 

n-1 refractivity 

p pressure 

q electronic charge 

R gas constant per unit mass 

R average R for the mixture 

R* gas constant per mole 

S fringe shift; also fringe shift at X = 4500A 

S' fringe shift at X' = 4030A 

T temperature 

v fluid velocity relative to the standing shock (lab) 

W atomic weight of argon atom 

W average atomic weight of the mixture 

x degree of ionization 

Greek Letters 

a polarizability 

7 ratio of specific heats 

A( ) change in ( ) 

5 shock stand-off distance 

9 argon ionization temperature 

X observing wave length; also X = 4500A 

X' 4030A 

p reduced mass of ion- electron pair 
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NOTATION (Continued) 



p mass density 

p mean density from region 3 to stagnation line 

X argon ionization potential 

u? frequency of observing light, corresponding to X 

ujp plasma frequency 

Subscripts 

1 region 1, initial downstream conditions 

2 region 2, conditions entering the standing shock 

3 region 3* conditions downstream of the standing shock 

s stagnation line conditions 

A atoms 

A ions 

e electrons 
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I . INTRODUCTION 



Treatment of high temperature shock layer flows has been of interest 
from the point of view of hypersonic and re-entry heat transfer. Fay 

r i*) 

and Riddell LlJ developed a correlation formula for convective stag- 
nation point heat transfer rate in a dissociating diatomic gas, in 
conjunction with an experimental study by Rose and Stark [2]. Several 
investigators have applied a modification of this formula to stagnation 
point heat transfer in partially ionized argon. Reilly [ 3 ]; operating 
in a regime where radiative heat transfer was considered negligible, 
measured convective heat transfer rates slightly higher than predicted 
by the Fay- Riddell equation applied to an ideal gas. Rutowski [4] also 
determined that convective heat transfer was only moderately increased 
by real gas effects in the partially ionized argon. Rutowski and Ber- 
shader [ 5 ] found that, for a shock Mach number of twelve, radiative and 
convective heat transfer rates were about equal, while the convective 
rate alone agreed with the Fay-Riddell ideal gas calculation. They were 
also able to correlate the measured radiative heat transfer to the cal- 
culated electron density, using the Kramers classical theory to relate 
radiation to the square of electron density. 

The importance of electrons in both radiative and convective ther- 
mal transport leads to interest in electron density measurement. The 
present study deals with a measurement of electron density near the 
stagnation line of a two-dimensional shock-heated argon flow around a 
cylinder-plate (Fig. l) . The purpose of the work was to make a quanti- 
tative comparison between measured results obtained by optical inter- 
ferometry and the predictions of equilibrium flow theory, using the 
Rankine-Hugoniot and Saha equations. A noble gas is desirable as the 
working fluid, since ionization can be studied without the complications 
of dissociation. In particular, argon was chosen because of its moderate 
ionization potential and relatively high molecular weight, which enabled 
significant ionization to be obtained behind the standing shock. 

'Numbers in brackets subsequently designate References collected 
at the end of this paper. 
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The technique of optical interferometry depends upon changes in 
index of refraction of the test gas, and has long been used to map 
densities in non- ionized flows. But Alpher and White [6] demonstrated 
that interferometry at optical wave lengths could be used to measure 
electron densities in a plasma. Since the refractivity of atoms is 
relatively insensitive to wave length, while the refractivity of elec- 
trons is highly dispersive, measurement of fringe shifts at two wave 
lengths enables the electron contribution to be separated. This yields 
a direct measurement of electron density, independent of the state of 
equilibrium of the plasma. Klein [ 7 ] applied this method to two-dimen- 
sional flow around a model in an electromagnet ically- driven shock tube, 
but the flow entering the standing shock was highly non-uniform. Wong 
[8] has utilized this two wave length technique to study the relaxation 
zone behind the moving shock in argon. His work, using the same combus- 
tion driven tube as the present study, reinforces the conclusion drawn 
from Figs. 2 and 3 that flow entering the standing shock is quite uni- 
form. 

The present study is part of a long range program at Stanford 
dealing with the dynamics and properties of shock-heated plasmas. In 
particular, the present problem has direct application to ionization 
phenomena in the close environment of hypersonic vehicles, such as 
communications blackout. 
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II. EQUILIBRIUM FLOW THEORY 



A. Introduction 

In order to estimate the stagnation region electron density, it 
is necessary to calculate the state properties of argon at the stagnation 
line. Beginning with the state of the gas in the expansion chamber prior 
to the arrival of the incident shock wave (Fig. 1, region l), the one- 
dimensional state calculations are made across the incident and standing 
shocks, assuming adiabatic flow and the attainment of local thermodynamic 
equilibrium after each shock. The equilibrium assumption is justified for 
region 2 by the work of Wong [8], as well as the interferograms in Figs. 

2 and 3 of the present work. The computation is completed by assuming 
that the gas flows isentropically and incompressibly from the standing 
shock to the stagnation line. 



B. Gas Properties Across Shock Waves 

Properties across both shock waves are calculated using the Rankine- 
Hugoniot equations for a partially ionized, monatomic gas in conjunction 
with the Saha equation. The validity of the Saha equation for the 
equilibrium argon plasma has been demonstrated by Wong [8] and Cole [12]. 
Results of a computer program obtained by Wong were used to obtain the 
state of the gas in region 2. The equations used across the standing 
shock follow. 



The degree of ionization, x 



N 



x = 



„ A ♦ . 



is defined by 



where N. , N and N are number densities of atoms, ions, and 
A ’ A + 

electrons in any region of flow. The Saha equation, applied to the 
atom- ion-electron gas, is 



x 



B^(T)php e X//kT + i 

2B + (t)(2tt |i) 3y/2 (kT) 5 / 2 
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where p and T are the pressure and temperature, hp is Planck 1 s 

constant, x is "the ionization potential, and p is the reduced mass 

of the ion-electron pair. In order to make use of the Saha equation, 

it is necessary to approximate B. and B , the partition functions 

A A + 

of the atom and ion. Since the neutral atom of argon has about 60 
excited states near the ionization potential, the partition function is 
approximately 

_ 162,500 

B a (T) = 1 + 60e T 



where l62,500°K. is the ionization temperature. Since the temperatures 
involved are less than 15>000°K. , it is evident that the second term can 
be neglected. The argon ion has a ground state of multiplicity four, and 
a low level excited state of multiplicity two. Neglecting higher levels, 
the partition function of the ion is 

2060 

B (T) = 4 + 2e T 
A 

Bond [9] estimates that neglect of these higher levels and the second 
ionization of the atom will not introduce appreciable error below 
15>000°K. Using these approximations, the Saha equation becomes 



Ph; 



x = 



182,100 

T 



- 1/2 



(2 V p) 3//2 (kT) 5 / 2 



4\2+e 



2060 \ 
T 



+ l 



The equation of state of an ionized monatomic gas in thermodynamic 
equilibrium is 



p = (l+x) pRT = pRT 



where R is the average gas constant per unit mass of the mixture , and 
the masses of the ion and atom are considered equal. 

The one-dimensional equations of conservation of mass, momentum, 
and energy for flow across a normal shock are unchanged by ionization: 



The velocities (v) are given in the laboratory system, i.e. a system 
fixed with respect to the standing shock. These equations may be 
combined to obtain 



h 3 • h 2 " I (^ - ^) ( p 3 ■ p a) ■ < II - 1) 

If excitation energy is neglected, the specific enthalpy is given by 
h = e. + e. + — , 

tip' 

where e^ is the internal energy per unit mass due to translation, and 
e^ is the specific ionization energy. They are given by 

e t = l lT = i=" T = il" (1+x)T = 1 R ( 1+x ) T > 



and 



e, 

1 



x(N + N )x N ke N k0 

— = = e ~y v ~T = xR0 

P P m a tN e /x) 



y 



where R* is the gas constant per mole, W is the atomic weight of 
argon, W is the average mixture molecular weight, 6 is the ioniza- 
tion temperature, and m^ is the mass of the argon atom. Finally, the 
specific enthalpy is 



h = | (l+x)RT + xR0 



Substituting this in Eq. (il-l), the Rankine-Hugoniot equation becomes 



(II-2) 



+ 1 



Q 



h + 4 - 

p 2 \ i+x 2 y t 2 



A third equation, obtained by applying the equation of state on each 
side of the shock wave, is 




A fourth independent relation is found by eliminating v^ from the 
continuity and momentum equations to obtain 

2 P 2 P 3 (p 3 /p 2 ) " 1 

v 2 " p 2 (pypg) " 1 

The densities are eliminated from the preceding equation by substitution 
from Eq. (II-2) and the equation of state. The result is 




(II-4) 



To obtain this equation in terms of Mach number, the speed of sound in 
region 2, a^ > is required. For a monatomic, ionized gas, a^ is 
expressed as 



a 2 S(1+X 2 )T 2 ' 

2 

Dividing Eq. (II-4) by a^ , and solving for the pressure ratio p^/p^ ’ 
the result is 



F = B \ 3 + 5M 2 + 5 



n M V + 32 (x 3~ x 2^ 0M 2 

1_M 2 15 (l+x 2 ) T 2 



1/2 ‘ 



where Mg is the Mach number entering the standing shock. This 
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equation, along with the Saha equation, and Eqs. (II-2) and (II-3), 




are valid when applied to equilibrium regions on each side of a thin 
normal shock. Since the thickness of the non-equilibrium region in 
Fig. 3 is about half of the shock stand-off distance, the assumption of 
one-dimensional flow through the shock may be questionable. 

C. Stagnation Line Gas Properties 

If the flow from region 3 to the stagnation line is assumed adia- 
batic and non-heat conducting, then the total enthalpy is constant along 
any equilibrium streamline. In particular, along the central stream- 
line, the enthalpy is given by 



The integrated one-dimensional conservation equations used across the 
shock are inappropriate as the flow becomes two-dimensional approaching 
the stagnation line. But Euler 1 s equation along the central stream 
tube can be written as 



where neglect of viscous stresses is justified by the relatively high 
densities, hence Reynolds numbers, involved. The assumptions of adia- 
batic, non-conducting, equilibrium, non- viscous flow imply the assump- 
tion of isentropic flow. Integrating the previous equation, we obtain 




vdv + — = 0 



P 




If the flow is also assumed to be incompressible, then p=p^=consto , 
and the result is simply Bernoulli’s equation: 
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The assumption of incompressible flow is hardly tenable for the hot, 
ionized gas, even at the low Mach number involved. But, since Rutowski 
[4] estimates the density change to be small, Bernoulli's equation will 
be used in these calculations. 

With stagnation enthalpy and pressure known, the previously derived 

equations of specific enthalpy and state, and the Saha equation are 

applied at the stagnation point to solve for p , x , and T 

s s s 

These equations are 



h = ^ (i+x )RT + x R0 
s 2 s s s 



= 



s ( 1+x ) RT 
s s 



and 



x - 
s 



p h e 

s p 



182, 100 

3 . T s 



(2 1 t n) 3//2 (kT ) 5,/2 4(2+e' 2060//Ts ) 



+ 1 



- 1/2 



Results of the iterative solution were used to calculate N from the 

e 

s 

relation 



N 



P x 
s s 



e m. 

s A 
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III. APPLICATION OF INTERFEROMETRY TO ELECTRON DENSITY MEASUREMENT 



A. The Mach-Zehnder Interferometer 

A schematic of the interferometer is shown in Fig. 9« Monochromatic 
light from the source is collimated before entering the half- silvered 
beam splitter S^ . Here the light is split into two beams, which travel 
via the fully silvered mirrors M^ or M^ to be recombined at the 
second beam splitter S^ . The light is focused by the lens L^ , and 
the camera is placed at the focal plane of the test section. 

If the mirrors were all perfectly parallel, and the optical path 
lengths of the two beams exactly equal, one would expect reinforcement 
of the in-phase beams, resulting in a uniformly bright image at the 
camera. Now, if one of the mirrors is rotated slightly, a uniformly 
varying difference in optical path is introduced, resulting in alter- 
nating reinforcement and destruction, i.e., a fringe pattern. If now a 
plasma is formed in the test section, the index of refraction, and hence 
the optical path length, of one beam is changed, resulting in a fringe 
shift given by 

0 A nL A (n-l)L 

s = *T = I — 

where A (n-l) is the change in refractivity, L is the test section 
width, and X is the monochromatic observing wave length. The com- 
pensating chamber is introduced to compensate for the windows of the 
test section of the shock tube. 

Since the primary quantity desired is electron density at the 
stagnation line, the problem is to relate the observed fringe shift 
(from ahead of the standing shock, 2, to the stagnation line, s ) to 
the change in electron density. Conditions behind the moving shock, 2 , 
are taken from calculations by Wong [8], which were verified experi- 
mentally by him. 
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B. Development of Fringe Shift Equations 



In order to relate the observed fringe shifts to electron density, 
the refractivities of the neutrals, ions, and electrons must be related 
to thermodynamic quantities. 

The refractivity of the neutrals is obtained from the polarizabi- 
lity of the atom, o:^ , by the equation 

2-n- a 

n A ' 1 = ^ a A N A = — (W = WA ' 

where = (2ir Q! /m^) is the Gladstone-Dale constant, which varies 

only slightly with wavelength. The atom number density is related to the 
degree of ionization by 



N A - ^ f 1 -*) • 

Therefore, the change in refractivity of the neutrals between any two 
regions is 



A (n A -l) = KjA [p(l-x)] 



Similarly, the ion refractivity is given by 



n + - 1 = 2tt a N 
A A 



» V ! 



a. 



A AAA 

The number density of the ions is 



a.m.N 
m A A A A + 



<V 

m A N „ + 



N = N = £5 
A + 6 m A 



The change in refractivity of the ions is 

V 

A ( n - 1) = K. — A (px) 

A + ^ a A 
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The refractivity of the electrons is given by Wong [8] to be 



( v d = - 



1 

2 



2 

2 



1 



where iOp , the plasma frequency, is given by 



bv q 2 N 



= 



m 



kir q^ 

m m. 
e A 



(px) 






w is the observing light frequency, and where the assumption has been 
made that the electron-atom collision frequency is much less than the 
observing frequency. This assumption is valid for conditions of the 
present study. The change in refractivity due to electrons is 



A (n e -l) = - 2 A (px) = - 6.77 X 10 8 X 2 A (px) 

2m m A 7rc 
e A 

where 1 is in centimeters, and A (px) signifies the difference in 
the product of density and degree of ionization between any two regions. 

Assuming that the refractivity changes for the various particles 
are independent, the total change in refractivity is 

A (n-l) = A (n -1) + A (n - l) + A ( n -l) 

A A 6 

Recalling that the total fringe shift across, for example, the standing 
shock is given by 



s = £ A (n-D 

this fringe shift is expressed in terms of thermodynamic quantities by 



S = L 



P 3 (i -x 3 ) -° 2 a-x 2 ) + — (p 3 x 3 -p 2 x 2 ) 



O 

-6.77 x 10 X(p 3 x 3 -p 2 x 2 



(HI-1) 
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The ratio of polarizabilities of the argon ion and atom is taken as 




and the test section width for the present experiment was L=5-0 cm. 
Inserting these values and combining terms, the total fringe shift for 
X = 4500 A is 



S = I.78 x 10 4 [p 3 -p 2 -8. 86^X3 - 
where = .1592 crn^/gm has been 

K^q 3 o = *1609 cm^/gm , and the fringe 

S’ = 1.993 x 10 4 [p 3 -p 2 -7.13(p 3 x 3 

Combining the last two equations, the 
the argon plasma is 



P 2 x 2 )] , (III- 2) 

used. Similarly, for X'=4030 A , 
shift is 

- p 2 x 2 )] . (in-3) 

electron density in region 3 of 



N = N + 4.33 x io 17 (S' - 1.1275) , (III-4) 

3 2 

where X = 4500 A , X' = 4030 A , and S and S' are the measured 

fringe shifts at these wave lengths. With N now determined, either 

e 3 

of the fringe shift equations, (III-2) or (III-3), can be used to solve 
for p 3 . Choosing the equation for X = 4500 A , the result is 



p = p + 5-63 x 10' 5 S + 5.86 x 10 " 22 ( N - N 

3 2 l e e r 



The degree of ionization can now be determined from 

N m. 
e A 

3 “ ~ p 7 ~ ■ 



The last three equations are sufficient to determine electron density 
and degree of ionization at any point in the plasma from observed fringe 
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shifts between a known region and the region of interest. The fringe 
shift equations (III-2) and (III-3) may be used to predict fringe shifts 
from calculated plasma conditions, as in Fig. 
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IV. DESCRIPTION OF EXPERIMENT 



The experiment was performed using the Stanford aerophysics 
laboratory shock tube. Figure 8 is a schematic of the experimental 
arrangement. The 25 foot long driven section is constructed of extruded 
aluminum, with a two inch square inside cross section. The combustion 
driver uses a stoichiometric mixture of hydrogen and oxygen diluted with 
helium, which is ignited by a battery of spark plugs. The driver and 

driven sections are separated by a scribed aluminum diaphragm. Pumping 

-5 

of the driven section, to about 10 mm. Hg before filling, is accom- 
plished through a station at each end of the section. The test section, 
located 20 feet from the diaphragm, supports the two-dimensional cylinder- 
plate model. A pair of optical windows permits viewing of flow about the 
model, and three barium titanate piezo-electric transducers just upstream 
of the model are used to measure the incident shock speed. The dump 
tank at the downstream end is separated from the driven section by a thin 
aluminum diaphragm which breaks on arrival of the incident shock. 

The Mach-Zehnder interferometer (Fig. 9 ), located at the test 
section, is on loan through the courtesy of Lockheed Missiles and Space 
Company. The light source was obtained by discharging l6 kilovolts 
from a 7*5 microfarad capacitor across magnesium electrodes spaced about 
l/4 inch apart. Although no direct measurement of the spark duration 
was made, oscilloscope pick-up of the capacitor "ringing" indicated 
that the spark was probably near maximum intensity during the entire 
exposure time of five microseconds. A blast shutter and pinhole were 
originally placed at the focal point of the light source in an attempt 
to reduce film exposure by plasma radiation. However, Fig. 10 shows 
that, even at relatively low Mach number, radiation reaching the film 
was intense enough to wash out completely the fringes through the 
5890A filter. 

Since it had become evident that the mechanical blast shutter was 
not fast enough for this application, a Kerr cell was obtained, again 
through the courtesy of Lockheed. This electronic shutter, manufactured 
by the Kappa Scientific Corporation of Sierra Madre, California, was 
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equipped with a pulse- forming network that permitted a very reproducable 
shutter-open time of five microseconds. For the present steady flow 
application, this exposure time proved satisfactory. Unfortunately for 
this application, it was designed for use in the ultraviolet, with 
acetronitrile as the working fluid. The Kerr cell was triggered to open 
for five microseconds during peak spark intensity, and the shutter-open 
period was always between 12 and 17 microseconds after passage of the 
incident shock by the model. 

The flow field was split by interference filters passing 4500 and 
5890 angstroms located just ahead of the Polaroid camera. Simultaneous 
measurements of near stagnation-line fringe shifts through the two fil- 
ters were sufficient to determine the electron density. However, 
variation in intensity of the dark fringes passed by the 589OA filter 
(see Figs. 2,3) indicated the presence of a harmonic. Further, compa- 
rison of observed fringe shifts with Eqs. (III-2) and (III-3) showed 
that the "5890" fringes actually corresponded to a wave length lower 
than 4500A . Examination of a spectrum obtained by Wong through this 
second order 589OA filter, also using a magnesium spark source, showed 
a third order pass band at 4030A . It was concluded that peaking of 
the ultraviolet Kerr cell at about 4000A , as well as the low intensity 
of the magnesium spectrum at 589OA , resulted in a fringe pattern domi- 
nated by k' = 4-030A . The existence of harmonics in the fringe pictures 
indicates that the presence of some 589OA might cause a small increase 
in effective wave length above the value of k' = 4030A used in the 
calculations . 

An attempt was made to obtain total stagnation line heat transfer 
data using a platinum thick- film calorimeter, as described by Rutowski 
[4]. Initial efforts indicated that this phase of the program would 
require more effort than originally anticipated. Because of practical 
limitations of time, it was decided to concentrate remaining efforts on 
the interferometric studies. 
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V. DISCUSSION OF RESULTS 



Usable fringe shift data were obtained for both X = 4500 A and 
X' = 4-030 A on four runs. These were made at incident shock Mach 
numbers of 11.45, 11.8, 12.0, and 12.8. The lowest Mach number was 
obtained with an initial pressure (p ) of 10 mm. Hg. , while the others 
corresponded to p^=5 mm. All tests were made in argon, with the l/2 
inch diameter cylinder- plate model installed in the test section. Graphs 
of fringe shift, density ratio, electron density, and degree of ioniza- 
tion (Figs. 4,5,6 and 7) compare predicted and measured quantities in 
region 3, just downstream of the standing shock. Figure 7 also shows 
measured stagnation line degree of ionization and the corresponding 
prediction of equilibrium theory. 

Attainment of equilibrium in region 2, the flow entering the 
standing shock, was found by Wong [8] to occur about six microseconds 
after the incident shock (lab time) for the lowest Mach number of the 
present study, while the prediction of Petschek and Byron [10] is 
somewhat less. Since the present interferograms were taken between 12 
and 17 microseconds after the incident shock, it is reasonable that the 
flow entering the standing shock was in equilibrium. The interferograms 
in Figs. 2 and 3 show relatively straight fringes in region 2, supporting 
the conclusion of equilibrium there. The compensating effects of inci- 
dent shock wave attenuation and radiative cooling on the flow in region 
2 were discussed by Rutowski and Bershader [5]. Although the magnitude 
of these effects has not been determined, their effect probably was 
minimized by taking all the interferograms at the same time relative to 
the incident shock. Total available test time was estimated to be about 
100 microseconds for the fastest run. 

The interferograms at the two lower Mach numbers (Fig. 2a) show a 
sharp upward (positive) fringe shift across the standing shock due to 
change in refractivity of the heavy particles [Eq. (ill-l)]. These 
shifts were found to correspond not to the predictions of equilibrium 
theory [Eq. (III-2)], but to the estimate obtained by assuming degree of 
ionization frozen at the value in region 2, and equilibration of the 
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translational modes only. Predicted quantities based on these assump- 
tions were obtained from the perfect gas shock tables of Keenan and 
Kaye [ll], for 7=5/3 , and the results are compared with measured 
values of fringe shift, density ratio, and degree of ionization in 
Figs. 4, 5, and 7° Since measurements at the stagnation line showed 
only a small increase in degree of ionization from that in region 3 
(Fig. 7), the flow field from the standing shock to the stagnation line 
is concluded to be definitely out of chemical equilibrium, and in fact 
nearly frozen. Measured values of electron density in region 3, which 
differed only moderately from those measured at the stagnation line, are 
shown in Fig. 6 to be much less than the corresponding equilibrium 
predictions. 

The interferograms at the two higher Mach numbers (Fig. 2b) show a 
radical change in the flow field. The shock stand-off distance is re- 
duced by more than a factor of two. The initial upward fringe shift 
again roughly corresponds to the frozen prediction, with the larger 
shift at the lower wave length, as would be expected from Eq. (ill-l), 
for the heavy particles. But a large shift in the opposite direction, 
i.e. due to electron production [Eq. (ill-l)] is also observed. Since 
this shift is proportional to wave length, it is greater for X = A500 A . 

The flow appears to approach local equilibrium at the line of 
maximum downward fringe shift, located about half way from the shock to 
the stagnation line. Since the fringe shift due to electrons from 
region 2 to 3 is somewhat less than predicted for equilibrium (Fig. A), 
the measured values of electron density obtained from Eq. (III-A) only 
approach the equilibrium estimate. It is possible that application of 
the one-dimensional conservation equations across a shock extending half 
way to the stagnation line is invalid, which could explain the lower 
measured equilibrium densities. A better theory would make use of a 
one-dimensional analysis across the frozen portion of the shock, and 
calculate the approach to equilibrium by a two-dimensional theory utili- 
zing a rate equation. However, the change in the fringe pattern (Fig. 3) 
perpendicular to the flow direction is small near the centerline of 
flow, indicating that the flow across the shock is nearly one-dimensional 
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in the region of interest. The flow from region 3 to the stagnation 

line causes an upward shifting of fringes, which for these conditions 

can only be caused by ion-electron recombination [Eq. (lll-2)j. This 

effect, which of course is not predicted by the adiabatic equilibrium 

analysis, is illustrated by a small decrease in degree of ionization 

from region 3 to the stagnation line at M =12.0 , and a larger decrease 

s 

at the higher Mach number (Fig* 7)* Although the Kerr cell effectively 
reduces the interferogram exposure from plasma radiation, other pictures, 
using only the blast shutter, show the existence of strong radiation in 
this Mach number range. It appears to be concentrated in the region 
behind the bow shock, and is probably largely due to radiative recombi- 
nation. It is concluded that the flow at the higher Mach numbers does 
approach equilibrium momentarily, but immediately begins cooling due to 
radiative losses. Radiative recombination could also explain the lower 
than expected equilibrium electron densities measured in region 3* The 
sharp change in stand-off distance over a very small Mach number range 
may be related to the attainment of equilibrium at the higher speeds. 
Table II summarizes these observations on the flow interf erograms . 

Although no rigorous theory has yet been developed for convective 
stagnation region heat transfer in a partially ionized monatomic gas, 
the Fay-Riddell correlation equation [l] has been applied to this situ- 
ation by Reilly [ 3 ] and Rutowski and Bershader [ 5 ]. Predicted heat 
transfer rates based on this equation depend on the value chosen for 
the Lewis number, but do not vary markedly from the ideal gas prediction 
based on the same equation. Both Reilly and Rutowski and Bershader 
found substantial agreement of measured convective heat transfer rates 
with estimates based on the equation of Fay and Riddell. The present 
study suggests that the flow field ahead of Reilly 1 s model was probably 
nearly frozen, since his study used a much smaller model and lower 
initial pressures than the present work, and involved Mach numbers of 
12 or less. A similar comparison of the present study with the experi- 
ments of Rutowski and Bershader, involving Mach numbers greater than 
11. 7 , initial pressures of 10 mm. Hg. , and a larger model, indicates 
that equilibrium was probably attained in the stagnation region during 
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their work. These considerations suggest that the state of equilibrium 
of the shock-heated argon plasma does not critically affect stagnation 
region convective heat transfer. This result is not entirely unexpected, 
since convective heat transfer is basically a, boundary layer phenomenon. 
In the first approximation, it should be relatively insensitive to shock 
layer rate processes, particularly for the relatively high densities 
involved in the present study (Table I). 

Rutowski and Bershader suggest a radiative heat transfer model 
based on an extension of Kramers classical theory for free- free 
(bremsstrahlung) radiation to include radiative recombination into the 
upper bound states. Assuming thermodynamic equilibrium they estimate 
radiative heat transfer rate to be proportional to the square of elec- 
tron density. Use of predicted electron densities in the present case 
would have resulted in grossly overestimated radiative heat transfer 
for the frozen tests, but a much smaller error at the higher Mach 
numbers, where equilibrium was attained. 

The following factors were considered as possible contributors to 
error of the experimental measurements: 

1. Error in determination of Mach number by the pressure 
transducer-oscilloscope system is estimated to be about 1$. 

2. Actual initial temperatures varied from T^=297°K to 

T^=304°K . Since al 1 theoretical calculations were based on 

T^=300°K , a small scatter in the (uncorrected) measured 
values would be expected. 

3. The work of Wong [8] indicates that the presence of impurities 
could cause a decrease in the relaxation time. However, this 
effect was predicted primarily for very low initial electron 
densities, where atom-atom collisions are predominant. Since 
all regions of interest in the present study involve relatively 
high electron densities, and the tube pressure before filling 
was about 10 mm. Hg. for all runs, the effect of impurities 
is not considered significant. 
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4. Fringe shifts could be read to an accuracy of about one-tenth 

of a fringe. An error in this measurement would introduce 
only a small discrepancy between predicted and measured fringe 
shifts at a single known wave length (Fig. 4, X=4500A) . 

Therefore, it is considered that qualitative observations from 
Fig. 4, i.e. separation of frozen and equilibrium regimes, are 
valid . 

5. However, the measured values of electron density, degree of 
ionization, and density ratio, which were all obtained from 
measured fringe shifts at two different wave lengths, are 
subject to more serious error. The proximity of the two ob- 
serving wave lengths (X=4500A, and X'=4030A) results in 
sensitivity of Eq. (III-4) to small errors in reading the 
fringe shift, especially at lower Mach numbers where the 
shift due to electrons is small. It will be recalled that 
original plans had called for using X'=5890A as the second 
wave length, which would have resulted in considerably better 
accuracy. There is an additional uncertainty in the actual 
value of the lower wave length (X'=4030A + ~100A , - OA) . 
This is derived from the possibility that presence of traces 
of 589 OA might result in an "effective" X' > 4030A . 
Elimination of all 589 OA radiation by procuring a suitable 
additional filter was precluded by time considerations. Use 
of a higher "effective wave length vice X'=4030A in the 
calculations would have resulted in a fairly uniform increase 
in the measured electron density for all runs. The estimated 
uncertainty in measured electron density due to these consi- 
derations is included in Fig. 6 . It should be noted that this 
interferometric method is capable of considerably greater 
accuracy than the present work obtained. However, the agree- 
ment of density ratios (Fig. 5) obtained by the two-wave 
length technique with the appropriate theory (frozen or equi- 
librium) corresponds to the conclusions drawn from the more 
accurate single-wave length fringe shift predictions (Fig. 4). 
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VI. CONCLUSIONS 



1. Flow entering the standing shock was in equilibrium during each 
test . 

2. Stagnation region equilibrium was not attained during the two 
lower Mach number tests. The degree of ionization remained 
essentially frozen at its value entering the standing shock. 

3. During the higher speed tests, local thermodynamic equilibrium was 
achieved ahead of the stagnation point* The thickness of the re- 
laxation region was about half of the shock stand-off distance. 

The degree of ionization after relaxation was less than the equi- 
librium prediction, and actually decreased from the end of the 
relaxation zone to the stagnation line* The lower than expected 
degree of ionization in both regions 3 and s is attributed to 
radiative recombination. 

k. Since radiative heat transfer has been related to the square of 

electron density, the present work indicates somewhat less radia- 
tive heat transfer than predicted by equilibrium theory* But the 
observed recombination provides some justification for considering 
radiative recombination a primary component of stagnation region 
thermal transport. 
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TABLE I 

CALCULATED STATE PROPERTIES OF EQUILIBRIUM ARGON 

-5 2 

Initial Conditions: = 5mm. Hg. , = 1.07x10 gm/cnr 

T = 300°K 



State 

Property 


Shock Mach Number, M 

s 




11 


12 


13 


P 2 (mm) 


795.5 


972 


1168 


P 2 (gm/cnr 5 ) 


4.94xio" 5 


5.54xio' 5 


6. 15x10" 5 


t 2 (°k) 


10,100 


10,880 


11,500 


X 2 


.0204 


.0305 


.0602 


N (cm -3 ) 

e 2 


1.52X10 16 


3 • 22X10 16 


1 /T 

5.59X10 


P 3 (mm) 


2534 


3642 


4990 


P 3 (gm/cm 3 ) 


12. 5X10" 5 


16.1x10' 5 


20. 2x10" 5 


T 3 ( ° K) 


12,065 


12,930 


13,820 


x 3 


.0705 


.1100 


.152 


N (cm -3 ) 

e 3 


1.33X10 17 


2.72X10 17 


4 . 62xio 17 


P s (mm) 


3090 


4340 


5448 


P 6 (gm/cm- 5 ) 


l4.8xio" 5 


18.8x10" 5 


21. 3X10" 5 


T (°K) 
s 7 


12,410 


13,240 


14,080 


X 

s 


.080 


.118 


.162 


N (cm -3 ) 

e 

s 


17 

1.79X10 1 


3 . 34xio 17 


17 

5.20x10 1 
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SUMMARY OF TEST OBSERVATIONS 
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X’=U030A 



X=4500A 



Fig. 2a. Flow Interferogram (M =11.^5) • 



X'=4030A 



X=^500A 





Fig. 2b. Flow Interferogram (M =11.8) . 

s 
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Fig. 3a. Flow Interferogram (M =12. 0) . 

s 




Fig. 3t>* Flow Interferogram (M =12.8) . 

s 
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FRINGE SHIFT, S , IN UNITS OF FRINGE SPACING 




Fig. k. Fringe Shifts. 
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A = 5 mm Hg 



Dp, 



10 mm Hg 



DENSITY RATIO ACROSS STANDING SHOCK 
vs 

INCIDENT SHOCK MACH NUMBER 

X = 4500A 
X' = 4030A 





11 



_J I I I L 

12 13 

M g (INCIDENT SHOCK MACH NUMBER) 

Fig. 5 . Density Ratio. 
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(cm iQ 1 ' 



0 



A = 5 ram Hg 



[D p^ = 10 mm Hg 



ELECTRON DENSITY IN REGION 3 
vs 

INCIDENT SHOCK MACH NUMBER 
X = 4500A 

X* = 4030A 



10 ' 



1 6 



11 



12 13 

Fig. 6. Electron Density. 
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M„ (INCIDENT SHOCK MACH NUMBER) 

Fig. 7 . Degree of Ionization. 
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DUAL FILTER 
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Fig. 9* Schematic of Optics. 





Fig. 10. Blast Shutter Interferogram. 
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